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Abstract 

We consider a membrane both weakly and strongly adhering to a ge- 
ometrically structured substrate. The interaction potential is assumed to 
be local, via the Deryagin approximation, and harmonic. Consequently, we 
can analytically describe a variety of different geometries: as well as ran- 
domly rough self-affine surfaces, smooth substrates interrupted by an iso- 
lated cylindrical pit, a single elongated trench or a periodic array of trenches 
are investigated. We present more general expressions for the adhesion en- 
ergy and membrane configuration in Fourier space and find that, compared 
to planar surfaces, the adhesion energy decreases. We also highlight the 
possibility of overshoots occurring in the membrane profile and look at its 
degree of penetration into surface indentations. 
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1 Introduction 



The statistical mechanics of membranes is an important branch of soft-condensed 
matter physics, not least because of its application to biological systems. Ex- 
amples of membranes that can be studied experimentally include liposomes or 
vesicles, microemulsions, lamellar liquid crystals, as well as biological cells, such 
as red blood cells [|1[|. Of great importance is a detailed understanding of the ad- 
hesion between two membranes or between a membrane and a substrate. This 
occurs ubiquitously in nature, with vesicle adhesion playing a dominant role in 
endo- and exocytosis ||, that is the communication of a cell with its immediate 
environment. Efficient drug delivery is dependent on the adhesion between a lipo- 
some and the plasma membrane of the target site while adhesion phenomena 
are also indispensable to biotechnology with, for example, biosensors being based 
on the binding of membranes to substrates. 

In this paper we choose to concentrate on the latter, that is the adhesion 
between a membrane and a solid substrate, in order to provide theoretical support 
for recent experiments aimed at creating new biotechnology. All of this research 
has involved the study of adhesion on materials which are not flat and chemically 
homogeneous but that have a deliberate chemical or geometrical patterning. As far 
as we are aware, there has been no theoretical work describing membrane adhesion 
on such structured surfaces. Previous studies have looked at vesicle adhesion on 
simple, planar, chemically homogeneous substrates f| and we hope now to 
redress this imbalance. 

First of all, we distinguish between weakly and strongly adhering vesicles or 
membranes. Weakly adhering giant vesicles have been studied using reflection 
interference contrast microscopy and have been found to lie at distances between 
300 A and 400 A from a planar substrate. Here, Van der Waals attractions to 
the substrate are mainly counterbalanced by a repulsive entropic Helfrich force. 
Strongly adhering or supported membranes J7J, however, sandwich a water or 
polymer film between them and the substrate and typically have a much smaller 
separation lying at distances between 10 A and 40 A. Hydration forces are now 
the dominant repulsive interaction. 

Supported membranes provide perhaps the most potential for biotechnological 
applications. They can be formed by the spreading of a bilayer over a substrate, 
vesicle fusion taking place at a substrate or by lipid monolayer transfer using a 
Langmuir-Blodgett technique These membranes are useful as they enable bio- 
functionalization of inorganic solids and provide a means to immobilize proteins, 
e.g. lipid coupled antigens and antibodies, with a well defined orientation and in 
non- denaturing conditions ||. This then offers a suitable environment to investi- 
gate protein-membrane coupling and protein-protein recognition processes. They 
also can be used to design phantom cells which allow study of the interplay between 
specific (lock and key) and universal forces during cell adhesion |]]. 

However, we wish to highlight the ability of supported membranes to act as 
biosensors. They can provide ultra-thin, highly electrically resistant layers on top 



2 



of a conducting substrate. If protein receptors are incorporated into these layers 
then electrical or optical detection of ligand binding to the receptors creates a 
biosensor. Alternatively, protein ion channels can be embedded in the membrane 
in which case their effect on its permeability endows the membrane with sensor-like 
properties. 

Proteins included in the membrane often have a dimension greater than its 
thickness (which is approximately 40A) and so can lift the membrane off the sub- 
strate. This is undesirable and can be prevented by creating water pockets on the 
surface of the substrate which then act as protein "docking" sites. Therefore, one 
is naturally led to consider membrane adhesion on geometrically structured sub- 
strates which have been indented in some way. This structure can be made by a 
variety of methods, for example, chemical etching of silicon wafers can lead, due to 
their particular crystal structure, to long "V" -shaped channels ||. Alternatively, 
one can use micelles made from diblock copolymers to coat any semiconductor 
wafer. If this is ion sputtered, then loading of the micelles before hand with an 
inorganic compound can provide large etching contrasts. Using this technique one 
can pattern a planar substrate with pits or islands on a nanometer scale [|K| . 

Motivated by such experiments, we discuss here a straightforward theoretical 
approach which allows the adhesion energy of a bound membrane to be calculated 
for a variety of different geometrical configurations of the substrate. The techniques 
used have been strongly influenced by recent work in wetting phenomena. There is 
now quite a substantial literature covering the effect of geometrical and chemical 
substrate structure on wetting films. A general form of Young's equation has been 



derived ]TI, I2| , the influence of disorder [|T^] and corrugation |TJ]] of the substrate 



on wetting transitions determined and the effect of chemical patterning [15| on 
the shapes of adhering droplets studied. However, we use as the inspiration for 
our approach, a systematic description of the configuration of wetting films on 
non-planar surfaces |16|, [T^, [18|] in which the full non-local form of the Van der 
Waals interaction can be accounted for. 

There are also many interesting applications involving the adhesion of mem- 
branes on chemically structured substrates. However, we will not discuss this 
further here but return to it again in an accompanying publication . 



To start, we describe the various molecular interactions included in our model 
and differentiate between strong and weak adhesion. In Sec. ^|, we consider simple, 
planar surfaces, which serve as basis for nearly all the following work. The most 
profitable analytic approach, the Deryagin approximation, is detailed in Sec. |. 
We consider (separately) a corrugated substrate (Sec. |Q| ) and one broken by a 
single (long) trench (Sec. fOP , a pit (Sec. |4.4|) and a periodic array of trenches (Sec. 
[4.5| ) . We calculate the equilibrium membrane profile and the adhesion energy. The 
variation of this energy with respect to the parameters characterizing the geometric 
structure is emphasized. In Sec. |4.6| , the adhesion energy for randomly rough or 
self-affine surfaces is given for two particular cases. Finally, we conclude with a 
short discussion. 
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2 The Free Energy 



To begin, we consider a membrane with an elastic modulus k and tension a inter- 
acting with a rough surface. For a free, almost flat membrane, which is infinite in 
extent, the bending energy can be described by an effective free energy functional 
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ds\ds2 y/g 



a + -K(2H) 2 



where (si,^) are coordinates in the membrane surface, g is the determinant of 
the metric, H = (c± + c 2 )/2 the mean curvature, and C\ and c 2 the two principal 
curvatures. As the Gaussian curvature is a total divergence, the Gauss-Bonnet 
theorem implies that it can be ignored for a membrane with fixed topology and we 
proceed to do just this here. Working in the Monge representation, let p = (x, y) 
be a two dimensional planar vector and the heights of the surface and membrane 
above some reference p-plane be z s (p) and h(p), respectively. The geometry and 
notation used is shown in Fig. [1]. 

For an experimental system, in which giant vesicles (rather than the infinite 
membrane considered here) are involved, the tension usually arises due to the 
conservation of the vesicle total surface area. However, an additional contribution 
develops, for a weakly adhering vesicle, as the lipids making up the membrane try 
and move from areas far from the substrate into the more energetically favorable 
region near the surface. However, for our case of an infinite membrane none of 
these intricacies arise and we can consider the tension as an external parameter. 

If the membrane thickness is S, the interaction between the surface and the 
membrane can be accounted for by a potential V(h, 5; z s ) and total membrane free 
energy is given as 



T[h] = J dpjyi + (V/i) 2 



a 



K 

2 



Vh 



l + (V/i) 2 , 



V(h,8;z s 



(2.2) 



21, E21. We will 



where we have used the Monge representation. 

The interaction can be broken down into several components 
always include an external pressure, a Van der Waals attraction and a repulsive 
force, whose nature will be different for strong and weak adhesion. In principle, 
other possible interactions, for example electrostatic forces, can be added for more 
refined calculations. The total potential, used in (|2.2j), is then 

(2.3) 



V(h, 5; z s ) = p ■ (h - z s ) + V vdw (h, S;z s ) + V iep (h; z s ) 



A different dependence on the surface height has been emphasized as, in general, 
the interaction potential is a functional of z s . Note that in ( |2.2j ), the potential V 
is integrated over the projected area, A = J dp, only because all non-local effects 
of the solid roughness and membrane fluctuations are, in principle, incorporated 
into the potential itself. 
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The existence of an external pressure p in ( |2.3| ) implies that there is a difference 
in pressure across the membrane. Recalling that the membrane is infinite in its 
lateral extent, the pressure couples linearly to the membrane height and could arise, 
for example, from the existence of macromolecules on one side of the membrane 
only. An additional contribution arises due to the finite membrane thickness but 
as we assume that this is fixed throughout the membrane, it is just a constant 
which we neglect. 

The Van der Waals attraction that the membrane feels towards the solid surface 
is given by 



W(h;z s )-W(h + S;z s ) 



(2.4) 



where W(h; z s ) is an interfacial-like Van der Waals potential. Equation ( pT4 ) arises 
due to the bilayer nature of the membrane and can be understood as a supposition 
of two wetting layers. One of these has height h + 5 while the other is of height h 
and has a Hamaker constant of opposite sign. For a wetting fluid film of thickness 
h(p) resting on a solid substrate, with a rough surface configuration given by 
z = z s (p), W(h; z s ) satisfies ]16| 



W(h;z s 



with a kernel interaction 



r rz s (p+p>) 

dz / dp' / dz w(p', z — z) 

h(p) J J-oo 



(2.5) 



A 

72 



2m+2 



"'( i' ) — I - I (2.6) 
n z \r J 

Following convention, we have denoted the Hamaker constant by A and will usually 
set the integer m = 2 to model non-retarded Van der Waals interactions. Equation 
( f2.5| ) comes from the sum over all possible pair interactions between the molecules 
making up two half spaces which are each capped by either the surface z = h(p) 
or z = z s (p). 

Notice that ( |2.5|) is a function of h(p) but a functional of z s (p). For a flat 
membrane on top of a smooth and planar solid surface, which will be discussed in 
the following section, this functional dependence on z s (p) can be ignored and then 
(for m = 2) 

w ^ = w.-ihh$ (2 - 7 » 

This is the standard result for the Van der Waals potential between two semi- 

apart 



22 



infinite bodies with planar surfaces held a distance h — z s 

To counterbalance this attractive force, a repulsive interaction is included and, 
as before, we define two different possible scenarios; the weak and strong adhesion 
regimes: 



2.1 Weakly adhering membranes 

In this case we consider a repulsive force which is entropic in origin. This arises 
due to the surface cutting off the region that can be sampled by the membrane 
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as it undergoes thermal fluctuations. A renormalization group description is the 
only general approach that can be used to describe such effects but, providing 
fluctuations do not become critical, a simple supposition ansatz can suffice |2l[| . 
In this, so-called weak fluctuation regime [p^j , an entropically based term is added 



to the membrane potential. Such a term was first described by Helfrich p4| for 
a membrane (with zero tension) in a stack of membranes and can be thought of 
as an ideal "gas" pressure. The presence of another membrane above and below 
introduces a new lengthscale describing the maximum wavelength at which any 
membrane within the stack can fluctuate without being strongly suppressed by 
its neighbors. Fluctuations with wavelengths below this lengthscale are essentially 
unhampered and so one can think of an ideal gas of uncorrelated membrane humps 
characterized by a single size. The pressure that this gas exerts is simply the 
Helfrich repulsion. 

When a membrane is under tension, the extent to which it can fluctuate is 
much reduced. Therefore, in some sense, the membrane "feels" the presence of a 
surface or another membrane less and there is a corresponding decrease in entropic 



repulsion. This effect was looked at by Helfrich and Servuss p3 |, who, within a 
mean-field theory, found that the entropic repulsion, between two membranes, de- 
cayed exponentially in h 2 in the tension-dominated regime. Such a result has also 



been confirmed by Evans [26] via a self-consistent mean- field method. However, 



by using a functional renormalization group approach ET[ to include thermal flue 



tuations, it has been argued that this decay, while remaining exponential, should 
be so in h only. This has been verified by Monte-Carlo simulation p7| . Analyti- 
cal forms of the entropic repulsion are only available for the tension and rigidity 
dominated regimes but ideally one would like to have a function which tends to 
these two expressions in the appropriate limits. While this function has been found 
numerically |27|, there does not yet exist a complete analytical expression for it. 



However, a simple argument by Seifert p] gives a form of the entropic interaction 
which takes the correct limit for a — > or h — > and has the dominant exponential 
decay in h also present in the tension regime. We choose, for simplicity, to adopt 
this form which is given by 

6T 2 / Q \ 2 

v ^ = JW~ K (shiSi) < 2 - 8 > 

where 

We denote the temperature by T and have set here and hereafter the Boltzmann 
constant ks to unity. The prefactor of (|2.8j) has been chosen to be consistent with 



the renormalization group description [21] as h — ► oo but we note that the exact 



value has no special importance in the model presented here. In the limit of zero 
tension (a/ka ~ a/T — > 0), ( |2.8|) tends to the well known result 3T 2 / (2n 2 K,h 2 ). For 
large Qh, the potential decays exponentially as exp(— 2Qh). In this limit, which 
can correspond to large a, the membrane should take on more interfacial-like 
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behavior and, indeed, the entropic repulsion for an interface in three dimensions 
is also expected to decay exponentially with height |23| . 

As 5 is constant we can ignore its influence on the fluctuation interaction; it is 
only the lower lipid leaflet which defines the volume between the membrane and 
the substrate and so determines Vfl uc . 

Equation ( |2.8| ) was initially assumed only to be valid for a membrane fluctuat- 
ing near a flat surface. However, a number of different methods have shown that 
the Helfrich term can also be used for a rough substrate if h is replaced by the 
local height of the membrane. This approach has been adopted by several authors 
considering spatially inhomogeneous scenarios [25, and agrees with scaling 



arguments based on exact solutions [[H|]. For temperatures nearing the unbinding 
temperature, where the thermal fluctuations of the membrane go critical |3"2]] , the 
validity of ( |2.8| ) can be questioned and so throughout this paper we will only con- 
sider temperature regimes far from the unbinding transition. Consequently, for a 
rough surface, we choose a fluctuation repulsion satisfying Vfl uc (/i; z s ) = V^ nc {h— z s ), 
with Vfl uc (/i) given by 



2.2 Strongly adhering membranes 

As the membrane now lies much closer to the substrate (a typical distance can be 
30A as compared with 340 A for the weakly adhering case), the repulsive interaction 
becomes dominated by hydration forces. While their exact origin is not entirely 
understood, see f22|, experimentally they decay exponentially j33f and have an 
Angstrom range, 



VhydW = bae 



-ah 



(2.10) 

where b is a dimensionless number and a an inverse length, typically a -1 ~ 2-3 A. 
For a rough substrate we again assume that we can replace h by the relative height 
and so use the local expression for the interaction potential, Vh y d(^ — z s ), in the 
free energy. Such an assumption seems reasonable due to the very short-range 
nature of the interaction in ( |2.10| ). 



2.3 The adhesion energy 

To summarize, the total free energy consists of 



a + -n{2Hf 



+p-(h-z s ) + V vdw (h,S;z s ) + V rep (h-z s ) (2.11) 



that is, a bending energy, made up of tension (coupled to the membrane area) and 
rigidity (coupled to its total mean curvature squared, H 2 ) terms, and an interaction 
potential V given by ( |2.3| ). The latter contains a contribution from an external 
pressure, the full non-local Van der Waals interaction and a local repulsion, V rep , 
which we consider either as a fluctuation, see ( |2.8| ), or a hydration interaction as 
in (]2T0|) . 



7 



One of the most relevant quantities in experiments is the membrane adhesion 
energy. Within our general mean-field approach, the optimal height of the mem- 
brane is that which minimizes fl2.ll ). The value of the free energy of the system 



when the membrane takes up this optimum configuration, J 7 ^, leads to a natural 
definition of the adhesion energy per unit area 

Here, A is the total area of the projected reference p-plane, A = J dp. 

Notice that in the above equation a tension dependent term is subtracted from 
the definition of the adhesion energy. Doing so conveniently shifts the origin in 
such a way that the adhesion energy for a completely flat membrane has no tension 
dependent contribution. This agrees with one's intuitive picture and is necessary 
as the membrane tension couples to the entire membrane area and not just to any 
excess area arising from a non-planar configuration. Consequently, for a membrane 
adhering to a flat wall and being itself flat, h(p) = ho for some constant ho, 
the adhesion energy is simply given as the negative of the interaction potential 
experienced by the membrane. For example, from ( J2 . 1 1| ) , [see also fl2.2p ], 

F min = Ao[o- + V(ho,5;0)] (2.13) 

and so ( [2.12| ) implies that 



U=-V(h ,5;0) (2.14) 

and the tension contribution vanishes. Hence, U is positive for all sufficiently 
attractive potentials, V. 



2.4 Rescaling of lengths and interactions 



Before we proceed and calculate the adhesion energy U for different types of cor- 
rugated and rough surfaces, it proves profitable to extract two natural lengthscales 
present in the problem. 

The first of these is provided by the ratio between the Hamaker constant, A, 
appearing in the Van der Waals potential (|2.5|), and the tension a p4|| , 



A y 



2trx 



(2.15) 



while the second describes the crossover between the tension and the rigidity dom- 
inated regimes 

£ = J± (2.16) 



Using these definitions ( |2.2|) becomes 

1 



l -m-jd P { 



1 + (Vh)< 



Vh 



1 + {Vh) 2 



+ -V(h,5; z s ) 
a 



(2.17) 



S 



Equation ( |2.5| ) gives 



-W(h;z s 



2a 2 r°° 



dz / dp' 



dz 



P 2 + (z~ z'f 



-(m+1) 



(2.18) 



O 7T Jh(p) 

while the fluctuation and hydration interactions can be written, respectively, as 



1 IT 

-v Ruc (h) = ^—^'- 2 

a Ann 



sinh- 2 (m) 



and 



-V hyd (h) =be 



-ah 



(2.19) 
(2.20) 



where f2 = {ixo /2T) l l 2 was defined earlier. 

In the next sections (|2.17 ) will be minimized and the adhesion energy U will 
be calculated for various corrugated and rough surfaces. However, first we discuss 
the case of a simple, flat surface. 



3 A planar surface: choice of physical 
parameters 

For a planar and homogeneous solid surface, mean-field theory predicts that the 
membrane also adopts a flat configuration, 



Hp) = h c 



(3.1) 



where we use the subscript zero to denote all planar quantities. Here, we set 
z s (p) = in order to ensure a zero average surface height. The Van der Waals 
term (for m = 2) also simplifies considerably 



a 



:Kdw(M;0) 



6 \h 2 (h + 6) 2 



(3.2) 



while the repulsion satisfies (|2.19|) and ( [2.20D for weak and strong adhesion, re- 
spectively. 

The membrane height ho obeys V = ^ = (balance of forces), that is 



V'(ho 



P+ 3 U§ (h + 5)* 



+ v: ep (h ) = o 



(3.3) 



with an adhesion energy, ( 2.12|) , of 



Un 



a 



V (h ,S) 



a 





a 2 


' 1 


i 


6 


M 




1 - 


a 2 


i 


1 


6 




(h Q +S) 2 \ 



+ j^sinh 2 (fi/io) weak adhesion 



+ be 



-aho 



(3.4) 
strong adhesion 
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where Vo(h,5) = V(h,5;0). Equation (|3.4j ) is useful as it serves as the starting 
point to which all our perturbation theories provide corrections. 

At this point in order to facilitate comparison with experimental systems we 
discuss the numerical values of our model parameters. A particular example is 
detailed in Table p], where T = 4.1 x 10~ 21 J at room temperature and kg was set 
to unity. For simplicity, the external pressure is set equal to zero throughout 

p = (3.5) 

Let us discuss now the choice of parameters for weak and strong adhering 
membranes, separately: 

(i) For weak adhesion, an effective value of the Hamaker constant is used, 
A = 8.67 x 10~ 22 J ~ 0.21 T. This value is quite small in order to approximately 
model the screening effect of ions in the solution surrounding the membrane |22[ . 
Setting a = 1.7 X 10~ 5 Jm~ 2 , the rescaling length, a ~ 28.5 A, and h is calculated 
to be ho — 11.85a ~ 338 A, in agreement with experiment ||. Such a value is 
reassuring as the expression, Q3.2| ) for the Van der Waals interaction, is only valid 
for h < 500 A before retardation effects begin to become important ||32|| . From 
Fig. [2], one can see that V (h , 5) ~ -1.1 x 10 _4 cr and so the adhesion energy 

U = -V (ho,S) 

~ 1.1 x 10~V = 1.87 x 10~ 9 Jm" 2 (3.6) 

which is positive as expected. 

(ii) If we now turn to strong adhesion, ion effects can be ignored and we use 
a larger value of the Hamaker constant, A = 2.6 x 10~ 21 Jm~ 2 ~ 0.63 T [0. This 
implies that a ~ 49.3 A and h ~ 0.61a ~ 30 A in agreement with measured values 



using specular reflection of neutrons The various parameters specifying the 
hydration force, see ( 2.10Q , are 



b = (0.93 Jm- 2 )/a ~ 5.47 x 10 4 ; a' 1 = 2.2 A (3.7) 



which are in accordance with those measured in [Q. This time, see Fig. §, 
V {h ,6) ~ -0.298(T and so 

U ~ 0.298a = 5.07 x 10~ 6 Jm" 2 (3.8) 



which is significantly greater than (3.6) 



Throughout the paper, we will keep to the particular values of the membrane 
and external parameters specified here and in Table [1]. It is sometimes convenient 
to express lengths in terms of the length a and energies (per unit area) in terms 
of the tension a, which is arbitrarily taken to have the same numerical value for 
weak and strong adhering membranes. We should say though that our model can 
offer only qualitative, or at best semi-quantitative comparison with experiment. 
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4 The Deryagin approximation 

One of the most useful approaches, which provides good opportunity for analytic 
progress, is the Deryagin approximation f36|. First of all, the full non-local Van 
der Waals potential, ( |2.4| ) and (|2.5| ), is replaced by a planar potential which is 
simply a function of the local relative height coordinate. Since the pressure term 
is always local and V rep already has this form, the total potential is written as 



V{h,6;z s ) ~ V (h-z s ,5) 

= p-(h- z s ) + V vdw (h - z s , 5; 0) + V Iep (h - z s ) 



(4.1) 



where Vo(h, 5) = V(h, 5; 0) as before. Notice that the above expression corresponds 
to replacing z s (p + p') in ( |2.5| ) with z s (p), i.e. removing the functional dependence 
of the potential on z s . The resulting free energy is then expanded to second-order 
in h — h and z s . Since the equilibrium position of the membrane is given by setting 
the variation to zero, the first-order term in h and z s vanishes, yielding 



a 
with 



T[h\^ JdpU + ^Vh) 2 + ^e^ 2 h) 2 + lv Q (h ,6) + ^(h-h -zA (4.2) 



8P_ 
dh- 



■V(h,5;0) 



—a 2 a 



h=ho 
1 



K (h + sy 



and 



3a 



[2 + cosh(2fi/i )] sinh- 4 (Q/i ) 
aa 2 be- aho 



weak adhesion 
strong adhesion 



(4.3) 



(4.4) 



Writing 5h(p) = h(p) — h , the Euler-Lagrange equation, giving the value of 
Sh which minimizes ( f4.2p , is relatively straightforward 



(kV 4 - aV 2 + v) 8h(p) = vz s (p) 



(4.5) 



To solve ( (4.5[ ), it is convenient to convert to Fourier space. Defining for any 
function f(p) 

/(q) = J dpf(p)e^ P (4.6) 



we find 



Sh(q) 



with the correlation lengths, having the usual definitions |21 



(4.7) 



(4- 
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For small g£ CT , ( |4.7| ) implies that the membrane follows the rough surface. However, 
as q^ a increases, both the tension and the rigidity act to dampen this effect. 
Using (|2.12| , |4.2|) and ( [4.7D , one can write down the adhesion energy 



a Aq J a 



dq 



A aJ 2(2vr) 2 



(l + q 2 £ + q'C) |^(q)| 2 + |5 s (q)| 2 - 2Sh(q)~z s (-q) 



I MO) 



(4.9) 



By defining the excess in the adhesion energy to be with respect to the planar case 
[recall that U = — V (h ,S)] 



AU = U-U 

= U + V (h ,5) 



we obtain 



AU 



1 



dq 



Q 2 C + « 



4<r2 



Fs(q) 



(4.10) 



(4.11) 



as higher order terms in q have not been 



2(27T^) 2 

which is only strictly valid up to 0(q A ) 
included in our starting equation. 

Equation ( [4.1 1|) has buried within it several assumptions which we would now 
like to make more explicit. It is important to realize that ( |4.7|) is not necessarily the 
general solution to ( |4.5| ) but just a particular solution. The general solution itself 



can be found by adding ( |4.7|) to the homogeneous solution, that is the solution of 
( (4.5[ ) when z s is set to zero. This then contains the four constants of integration 



required to satisfy any boundary conditions. Therefore , ( 4.11 ) is only correct when 
the homogeneous solution of ([4.5| ) is identically zero, which is true, for example, 
at a sinusoidally corrugated substrate (see Sec. fO])- Otherwise, while ( [1.1 1| ) is 
certainly included in U, it is not the whole picture and in this case it may well be 
better to work entirely in real space. 

For the particular experimental system specified in Sec. |3], one can calculate the 
various correlation lengths. These have been included in Table [IJ for completeness. 

The approximations used in this section have, as their basis, essentially a per- 
turbation theory (assuming h — ho — z s is small) around the planar value of the 
adhesion energy, Uq. Generally, one can only believe such an approach if the cor- 
rection term, AU, is much smaller than the result it is trying to improve upon. 
Consequently, our findings are strictly only valid when 



AU 



U 



< 1 



(4.12) 



which in fact limits the roughness of the substrates we can consider. 
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4.1 Sinusoidally corrugated surface 

The simplest case to look at is a sinusoidally corrugated surface 

z s (p) = csin(gx) (4-13) 
with an amplitude c and period 2n/q. Solving (|4.5|), one can show that 



r, / x csin(gx) . A 1 A . 

"V> - rrmrk (4 ' l4) 



Similarly, AU in ( |4.ll| ) is 

AU c 2 / q 2 £ + qV K 



4e vi+? 2 e+? 4 e 4 



(4.15) 



which is valid for small eg. 

In Fig. |3|, Z7/Z7o (using the parameters of Table |l|) is plotted as a function of 
the rescaled wavenumber aq for both weak and strong adhesion. From (|4.14 ), the 



average height of the membrane is unchanged from the planar result h . However, 
as q increases and the surface becomes corrugated with shorter wavelengths, the 
adhesion energy decreases (relative to the planar case) due to the extra bending 
energy cost as the membrane tries to follow these surface configurations. For large 
q, U flattens out; any point on the membrane is almost equidistant from a crest or 
trough on the substrate surface and so there is little energetic benefit in mimicking 
them. From ( |4.l5[ ), we have, in the limit q 2> max^ 1 ,^ 1 ) 



AU ~ -— (4.16) 

which is in accordance with U/Uq ~ 0.54 and U/Uq ~ 0.51, see Fig. |3], for weak 
and strong adhesion, respectively. For small q <C min(^~ 1 , i^" 1 ), the membrane 
can always follow these long wavelength perturbations, 

AU ~ -^p (4.17) 

and U remains q dependent. 

It is worth pointing out the very different x-axis scales and choices of c in Fig. ^. 
For weak adhesion c = 5a, that is c ~ 143 A, the adhesion energy "bottoms out" 
for q ~ 0.03/a. For strong adhesion c is much smaller, c ~ 8 A and much larger 
values of q are needed before the structure of the substrate is effectively smeared 
away. The values of c in both cases were chosen to ensure that ( |4.12| ) remains 
valid and one can see that U/Uq is never much smaller than 0.5. Our perturbation 
theory seems to work for amplitudes, c, which do not grow significantly bigger than 
around one third of the height taken by the membrane above a planar substrate, 
i.e. ho/3. 
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By noting that £ 2 = ( |4.15| ) can be written as 

AU _c* / q 2 



(4.18) 



and so one can also look at the effect of the elastic modulus k on the adhesion 
energy by plotting U/U against the crossover length £/a while keeping £ CT constant, 
see Fig. [|. As £ increases, k becomes greater than the strength of the Van der Waals 
interaction, the Hamaker constant, A. This increased bending energy (coupled 
with the fact that the membrane average position is fixed at ho) leads to a raise 
in magnitude of AU and so the adhesion energy decreases relative to Uq. 



4.2 The profile equation for piecewise constant surfaces 

In some cases, and in particular if z s is piece-wise continuous, it is more convenient 
to work directly in real space (as opposed to Fourier space). The task is then to 
solve ( P~5| ) and use ( |4.2| ) to find U. Given a piecewise constant surface, ( |4.5| ) 
becomes 

(/tV 4 - aV 2 + v) (5h(p) -z B ) = (4.19) 

for a constant z s (but which changes value discontinuously for different ranges of 
P)- 

The operator in ( [4.19 ) can be factorized 

kV* - aV 2 + v = k (V 2 - ri 2 + ) (V 2 - r? 2 ) (4.20) 
where rj± are complex in general and satisfy 

l±v^4^" 1 



v± = r 1 

with 



(4.21) 



X={™/a 2 f = j- (4.22) 

SK 

a dimensionless ratio of lengths. For our chosen set of parameters (Table HD, 
X — 0.302 and x — 9-44 for weak and strong adhesion, respectively. 

Therefore, to solve (|4.19|) which is a 4 th order differential equation we can 
consider, equivalently, two coupled 2 nd order differential equations 

(V 2 - V l)(5h(p)-z s ) = u (p) (4.23) 
(V 2 -vl)u (p) = (4.24) 

for the function Sh(p) and a second function uo(p). However, because of the 
underlying linearity of the 2 nd order differential equations, the general solution is 
simply 

5h(p) = c u (p) + cmi(p) + z s (4.25) 
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for constant Co, C\ and z s , and where u\(p) satisfies 



(V 2 -7] 2 _) Ul (p) = (4.26) 

Thus, we have reduced the 4 th order differential equation to two familiar 2 nd order 
Helmholtz equations. This approach is most useful for systems with cylindrical 
or spherical symmetry. We note that for those with one- dimensional symmetry, 
z s (p) = z s (x), it is straightforward to solve the 4 th order differential equation ( |4.19|) 
directly. The general solution is (for constant c ,Ci,c 2 and c 3 ) 

6h(x) = c e^ x + Cl e-" +I + c 2 e^ x + c 3 e~^ x + z s (4.27) 

which we will use in the following sections. 



4.3 An isolated trench 

We next consider a single, infinitely long trench parallel to the y axis and of width 
d and depth Xd. The height profile of the surface takes either of two fixed values: 
—Xd for |x| < d/2, and zero otherwise (see Fig. [|). This can be also written as 



z s {p) = z s (x) = -Xd{9(x + d/2) - 6(x - d/2)} 



(4.28) 



where 9(x) is the Heaviside step function, being only non-zero (and equal to unity) 
for positive x. The one-dimensional solution of Q4.19| ) is given by ( f4 . 2 7| ) with z s 
taking either of the fixed values —Xd for \x\ < d/2 or zero otherwise. The constants 
of integration can be determined from the following boundary conditions 



5h(-x) 
5h(±x) 



5h(x) 




cLS OC 



oo 



(4.29) 
(4.30) 



and by imposing continuity (up to the third derivative in x) at the edges of the 
trench, x = ±d/2. The final solution is (for positive x > 0) 



-5h(x) 



-1 + lk-e~— 

1 7 



cosh(r] + x) 
cosh(?7_x) 



— hk- sinh ( ) e v+ 



— \k + sinh 



2 



-r\-X 



for < x < d/2 



for x > d/2 



(4.31) 



and 



k± = 1 ± ; (4.32) 

In Fig. |5], the membrane profile is sketched for a supported membrane (strong 
adhesion case) with the same parameters as those given in Table 0. The trench is 
specified by d = 3a, i.e. 148 A, and A = 0.2. The membrane can be seen to follow 
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the contour of the surface, as expected from the general prediction ( |4.7| ), without 
developing a similar discontinuity to that occurring at the boundary of the trench. 
In addition, it can be seen that an overshoot is present, with the membrane having 



a height greater than /i (f° r an extensive discussion of overshoots see [30]). From 
([131]), 5h = h-h o can be shown to be negative for all x > d/2, providing t]± are 
real, and so the overshoot can only arise if r)± go complex. That is, from ( |4.21|) , 
when 1 — 4% 4 < or 

2£ > £ CT necessary condition for overshoot (4.33) 

i.e. in the rigidity dominated regime (4k, > cr 2 /v). 

Equation ( 4.33|) can be understood by realizing that the rigidity term in the 



free energy (|4.2| ) prevents the membrane from turning any sharp corners (see Ap- 
pendix). For a pure interface (k = and o ^ 0), ( |4.33| ) cannot be satisfied and, 
within our Deryagin approximation, no overshoot ever occurs. 

It is also interesting to look at the extent to which the membrane penetrates 
into the trench. A natural measure of this quantity is the membrane height in the 
trench center: 

5h(0) = \d (-1 + -A;_e-"+ d + -k+e^-A (4.34) 

In Fig. |] this quantity is plotted against d/a, the trench width [here, A simply 
scales out - see ( |4.34| )1. For larger d, the membrane is able to enter the trench 



more easily as there is less bending energy cost in the smoother configuration 
required to do so. For small d/a, Sh(0) — > and h(0) tends towards the planar 
value of the average membrane height as the trench becomes increasingly more 
narrow. In the opposite limit (not yet visible in Fig. H), d/a 3> 1, (|4.34|) implies 
that Sh(0)/a — ► —Xd/a, though here, we are pushing our perturbation theory 
beyond its region of validity. 

The adhesion energy ( |2.12| ) is defined as an energy per unit of projected area. In 
order to obtain a finite contribution in the case of an isolated surface perturbation 
like a trench, we need to consider one embedded in a flat surface of finite lateral 
extent L. Hence, using a local definition, involving a cut-off L, we can write that 
the change in U, AU (up to second-order in Sh) as 

AU = --V [ L dx \a(5h'f + K(5h") 2 + v(5h- z s ) 2 } (4.35) 
2L Jo 1 ' 



From ( 4.27|) , we can see, that for real rj, the cut-off should be proportional to rj + l 



due to the exponential decay. However, for \ > l/v2 both 77+ and 77_ are no 
longer real and it is straightforward to show that 



1 



v+ = v- = ^r 1 (V2x 2 + 1 + *vV - l ) ( 4 - 36 ) 

where the asterix denotes complex conjugation. In this case, (|4.27|) implies that for 
x > d/2 the profile has an exponential term, exp ( — — x J , and so a natural 
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choice for L is 



8 x max 



'd 



(4.37) 



where the prefactor was chosen by examining some numerical solutions for the 
profile. If r/ + is complex then ( |4.36| ) implies that l/$t{rj + } = —f=L=. Using this 



V^~ 2 

definition L ~ 22.3a for the system specified in Table |], which is a sensible value 
(see Fig. ||). 

Using (|4.31|) and (|4.35|) , one can show that 



At; 



a 



\ 2 d{ k 



32L 



XT J ^ A +) + 




/(A_,A_) 



2£ 



+ MM M-/(A+,A-) 



(4.38) 



where A± = rj±d is introduced, 



I(u,v) = uv(u + v)(e u - l)(e v - 1) e - (u+u) (H#) 



x 



>+■»)( 



u(c v -l)+v(c u -1) 
ue u (e v -l)+ve v (e u -l) 



(4.39) 



and k± are defined in (|4.32|) . The function I(u,v) can be seen by inspection to 
be positive, given that 2L > d (from Q437|) ), implying that AU is negative. Any 
non-planar membrane configuration (and such a configuration must be adopted by 
a membrane adhering to a rough substrate) will give an additional (and dominant) 
bending energy cost in the definition (|2.12j ) and so decrease U. 

Equation ( 4.38Q is illustrated in Fig. [7[ The adhesion energy decreases for 
larger d as the membrane adopts a more and more non-planar configuration. For 
small d/a, one can show that fl4.38| ) becomes 



AU 



4£ 2 L 




(4.40) 



and vanishes as d — >• in agreement with Fig. |^. For larger values of d/a, Fig. [7| 
would seem to imply that AU continues to grow in magnitude. This, of course, is 
false and is an unfortunate artifact of going beyond the valid limit of our pertur- 
bation theory. Exact numerical solutions for a similar scenario will be discussed 
in a companion paper [HJ and do not exhibit such unphysical behavior. 

Again, we wish to emphasis that the parameters specifying the substrate ge- 
ometry (A and d in this case) were chosen so that the constraint ([4.12 ) was obeyed. 



4.4 An isolated pit 

We next consider a single cylindrical pit of radius r and depth Ar. Equation ( 4.19|) 
is now only a function of p = \/x' 2 + y' 2 and is solved by Bessel functions, as is 
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usual for systems with cylindrical symmetry, 



Sh(p) = c l {pv-) + ciK {pr]-) + c 2 Io{pv+) + c 3 K (prj + ) + z s (4.41) 

Here, z s = — Xr for p < r and is zero elsewhere. Requiring a finite solution at 
p = and a vanishing one as p — > oo, implies that 

§h r p \ = [ c oh(pV-) + c 2 Io(pV+) - Ar for p < r 
\ CiK Q (pri^) + c 3 K (pr] + ) for p > r 

The constants of integration can be found by imposing continuity at p = r. 

The cylindrical symmetry results in little qualitative changes from the previous 
section (see Fig. [8]). Though, our perturbation theory is now acceptable up to a 
large hole of radius r = 350a ~ 0.997 pm due to the much smaller non-planar area 
of the substrate. A quick glance at Table [I] shows that in this case rj± are real and 
from ( 4.21 ) and ( |4.33| ) no overshoot can then occur ||3~8| . This can also be seen in 



Fig. [8[ A local adhesion energy can be defined similarly to ( |4.35| ) with the cut-off 
L obeying 

L = 8xm -(''-*{b) < 4 - 43) 

Unfortunately, for this case AU is a very complicated expression and it is not 
possible to proceed analytically. Therefore, similar to Sec. fO| AU is calculated, 
via ( [4.35| ), with the solution ( [4.42j ) but then the integrals are evaluated numerically. 
One finds that U decreases with r in a similar manner to Fig. ^. 

4.5 Periodic array of trenches 

The next, more complicated, scenario to consider is a periodic array of one- 
dimensional trenches. These are of depth Xd and width d in the x direction, 
and infinitely long in the y direction. We let the length of the repeating unit 
making up the surface be pd (see Fig. |9|), then 

■y oo 

z s (p) = -Xd- Xd [0(x-npd + d/2)-9(x-npd-d/2)} (4.44) 

^ n=— oo 

summing over integer n. 

By ensuring that the solution is locally symmetric about each trench center, it 
is possible to write out, using ( |4.27| ), 5h explicitly over one period 

• for < x < d/2 

5h(x) = c cosh(^ + x) + c\ cosh(^_x) — — (4.45) 

• for d/2 < x < pd — d/2 

5h(x) = c 2 e v+x + c 3 e- v+x + c 4 e v ~ x + cge"^ + y (4.46) 
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for \id — d/2 < x < \xd 



5h(x) = Co cosh (77+ (x — fidfj + c% cosh^?7_(x — fid) 



\d 



(4.47) 



As the configuration taken up by the membrane is symmetric around x = fid/2, 

5h(x) = 6h(/id - x) (4.48) 

one can show that 



c 3 = c 2 e 
C5 = c^e 



r\-iid 



(4.49) 



and so reduce the number of unknowns to four. These can be found by imposing 
continuity of 8h and its derivatives at x = d/2. 

The change in adhesion energy (to second-order in Sh) is 



AU 



d/2 



dx 



[id Jo 

1 



a(5h') 2 + K(5h") 2 + v(5h + Ad/2) 2 
~a(5h') 2 + n(5h") 2 + v{5h - Ad/2) 2 



"\2 



: / dx 

fid J d/2 



and after some algebra, can be found to satisfy 

AU A 2 



(4.50) 



a 



2(JL 



/(A + ,A_) + /(A_,A + ) 



(4.51) 



where we define 



I(u,v) 



(e u - l)(e^" 1 ) u - 1) 



uv 



QflU _ X 



{u 2 — V 2 ) 2 



(4.52) 



recalling that K± = rj±d. Again, it is possible to show that AU is always negative. 
Note, that the dependence of the adhesion energy on the depth of the trenches, 
via A, is, from ( |4.51| ), simply parabolic (providing ( |4.12| ) holds). 

For a surface roughness of this form, an interesting consequence of ( |4.51| ) is 
that the excess adhesion energy, AU, has a minimum as a function of fi at /1 = fi*. 
This is illustrated for a strongly adhering membrane in Fig. [l(| If we use e = d/£ 
as an expansion parameter, (for Fig. |TD|, e ~ 0.11), then perturbation theory gives 



AU 

AV ~ 27^" 
which implies 



.4 2 



X e + 



0* - ^ 

1440^ 



(^ + 2/i-2) X 8 e 6 + 0(e 



2-^-e 4 
30 



O e 



(4.53) 



(4.54) 
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This can be understood by looking at the average height of the surface, 



dx z s (x) 




(4.55) 



For 1 < jj, < 2, as fi increases, (|4.55|) evaluated near fi = 1, implies that the 
surface can be visualized as an array of spikes perturbing an "initial" planar state 
with (z s ) = —Ad/2. For \i = 1 the substrate is flat and consequently the general 
solution, ( |4 . 2 7| ) , implies that the membrane is so also, with Sh = —Xd/2 and 
U = Uq. As the thickness of the spikes increases, the membrane responds and is 
forced to bend more and more. This costs energy and AU increases in magnitude. 
One can show that 



AU 



a 



Mr 



^-i) + o((^-i) 2 ) 



(4.56) 



for // ~ 1. However, for /i > 2, the surface (and hence the membrane) starts to 
tend towards another, "final", planar state, that of (z s ) = Xd/2 which occurs when 
/i = oo. As Sh of the membrane, at this extreme point, also satisfies Sh = Xd/2, 
\AU\ decreases with growing /i. Looking at the two extreme values, /i = 1 and 
fi = oo, one can see, from ( 4.51|) and ( |4.52|) , that AU vanishes as expected. 

This behavior becomes more transparent if the variance in the height of the 
surface is considered 



<A* S 2 > 



<(^s-^» 2 > 

iff-*-* 



(4.57) 



with the average defined in ( (4.55| ). Equation ( f4.57| ) vanishes for /i = 1 and /i = oo, 
and has a maximum at fi = 2. Therefore, we find that the excess adhesion energy 
is greatest (in magnitude) for that value of [i at which the surface is the most 
rough, i.e. \i = 2. In fact, if U is plotted against (Az^), it is just a monotonically 
decreasing function. From (|4.53| ) and ( f4.57|) , we have 



AU 



a 



X 4 (Az s 2 ) 2 
2d 2 



+ 0(e«) 



(4.58) 



and, at this level of approximation, there is a simple, linear relationship. The 
average height of the membrane itself can also be calculated 



(h) 



h + {Sh) 



h + 



Xd 



1 - 



(4.59) 



and /j, = 2 again arises as a significant value. 
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It is also interesting to examine the extent to which the membrane penetrates 
each of the trenches. Concentrating on Sh(0), we find 

Sh(0) a2a2 / J(A + ) - I(A_) 1 1 

-JdT = A+A - 1 Al-Al | " 2 (4 ' 60) 

where 

J( u ) = -—- ■ ! 1. (4.61) 

w m 2 e^ u - 1 

and A± = 77-1-d. Hence, at the level of the Deryagin approximation, there is a linear 

dependence on the trench depth, Ad. Looking at ho + Sh(0) — (z s ), from (|4.55|) , one 

can see (as I(u) is a monotonically increasing function) that the membrane always 

penetrates further into the trenches as their depth grows. It experiences a more 

attractive potential due to the greater surface area of the substrate. However, one 

must be aware that increasing A can quickly cause ( |4.12j ) to be broken. 



Finally, in Fig. |TT| we plot 5h(0) against the periodicity parameter, /i. It is, 
perhaps, at first sight surprising to see that for large /z, Sh(0) does not vanish. 
However, this is the wrong function with which to take this limit, as one always 
has at least one trench in the system and is always looking at the depth in this 
trench. If instead 5h(fid/2) is considered, this tends quite rapidly to Xd/2. For 
example, with the parameters given in Table |1] and Fig. [11], 5h(d) ~ 0.120a when 
H = 25. Hence, the membrane returns to its height above a planar substrate. 
However, one can see from Fig. [TT] that for « ~ 8, the trenches no longer have 
a significant effect and 5h(0) plateaus; the membrane effectively experiences a 
potential generated by a flat surface. 



4.6 Self-affine surfaces 



In this penultimate section, we describe surfaces whose roughness is random in 
character. This is normally used to model the "natural" roughness of a substrate 
which arises due to chemical impurities or crystal defects. We adopt the self-affine 
approximation [39| and assume that the substrate appears relatively smooth on 
large lengthscales. Typical fluctuations of length L grow only as where (3 is 
the roughness exponent < /3 < 1, which is directly analogous to the thermal 
roughness of an interface [ 23| . However, such models become divergent as L — > 00 



and we prefer to adopt a reduced self-affine scaling law |5J where a long lengthscale 
cut-off has been explicitly included 



Imp +p) - z s {p 



2 ) = 2 7 2 l-e 



(4.62) 



,2/3 



or equivalent ly 

(z s (0)z s (p)> = 7 V^ (4.63) 

Here 7 is the saturated root mean square roughness at long lengthscales, and ( is 
the crossover length from self-affine to saturated behavior. The wetting properties 
of such surfaces (k = 0) have already been quite extensively studied |13, [IH] . 
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To be able to calculate the membrane adhesion energy analytically, we special- 
ize to two cases: 



4.6.1 Roughness (3 = 1 

Equation ( 4.63|) can be Fourier transformed to give 

2 2 

(5 s (p)S s (q)> = 4tt 3 7 2 CV^(p + q) 
which can then be used in ( |4.11| ) to calculate AU. We find that 



AU 



a 



1 {~fX 



2\ 2 



where 



I(u) 
A± 



/(A+) - /(A. 
Ai - A 2 _ 



u 



and 



r(o,«) 

is an incomplete Gamma function. 

4.6.2 Roughness (3 = \ 

In this case we have 

(z s (p)z s (q)) = 



OO g 

dz— 



2vr 3 7 2 C 2 
(1 + CV) 



and 



with 



and 



AU 



a 



1 / 7 ^ 2 



1 + 



1 1 , 

I(U) = — + — lOE 



IT 



2u 3 



j5(p + q) 



/(9 + )-/(a 
e 2 - e 2 _ 



1 -u 
l + u 



©I 



(4.64) 
(4.65) 

(4.66) 
(4.67) 

(4.68) 



(4.69) 



(4.70) 



(4.71) 



1 - (C^ ± ) 2 (4.72) 

For 7 = 4a, equations ( |4.65| ) and ( |4.70| ) are illustrated in Fig. |I2]. As the 
crossover length, (, grows, the disorder becomes more and more correlated and 
the substrate becomes smoother on larger lengthscales. One can see that in both 
cases, U is tending towards the planar value. For smaller f3, this plateau is reached 
more quickly in accordance with ( [4.63[ ) . As ( decreases and the substrate becomes 
rougher, the adhesion energy deviates further and further away from the planar 
value. 
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4.7 Tensionless membranes 



So far we have only considered membranes which have a finite tension, a, and now, 
for the sake of completeness, we will briefly illustrate how our results change in 
the limit of o — * 0. For vanishing membrane tension, the definitions ( [2.16| , 4.21j) 
and (|4.22p imply that 



for cr = 



(4.73) 



Equation Q4.27| ) then indicates that the profile will have an oscillatory though 
exponentially damped form. 

To see how the adhesion energy changes, we can look at, for example, the 
adhesion energy of Sec. [4.5| , equation ( |4.51| ) which is valid for a surface periodically 
patterned with trenches. Multiplying ( [4.5 1|) by a, one is naturally led to consider 



(4.74) 



in the limit of zero tension. A little thought then gives 



AU 



(d_ 



with 



I(u,v) 



e « _ i)( e 0*-i)« - 1 

QflU _ I 



/(A + ,A_) + /(A_,A + ) 



uv 



(4.75) 



(4.76) 



{u 2 — v 2 ) 2 

and A± = dr]±, with the rfs obeying ( |4.73| ). While it is not as obvious as before, 
AU is again negative. 



5 Discussion 

In this paper, we have looked at the effect of geometric surface structure on the 
adhesion properties of a membrane and a solid substrate. The membrane interacts 
with the substrate via a Van der Waals potential and experiences either a Helfrich- 
like entropic repulsion, if it is weakly adhering, or hydration forces, if it is strongly 
adhering or supported. We find that the most analytically successful approach is 
that of Deryagin in which the entire interaction potential is assumed to be a local 
function of height. By making a further harmonic approximation quite a large 
number of surface geometries can be solved for. In general, the membrane tries to 
follow the surface contour of the substrate but for short wavelength undulations is 
damped due to its tension and rigidity. Two lengthscales, £ CT and £ K , analogous to 
the "healing" length |l6j for interfaces, emerge and control this damping. 

Within the Deryagin description, we find that the membrane can overshoot 
(similar to the behavior in a laser trap ||29|| ) near a disturbance in the planar 
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structure of the surface. Overshooting (i.e having a height greater than that it 
would have had for a completely flat surface) is seen to be a consequence of the 
eigenvalues of the 4 th order operator, making up the Euler-Lagrange equation, 
going complex. This occurs for 2£ > £ CT and consequently is directly due to the 
membrane having a finite (non-zero) rigidity. 

In general (but still at the Deryagin level), increasing the roughness of the 
surface reduces the membrane adhesion's energy. Having a non-planar substrate 
immediately leads to a competition between the bending and potential energy 
contributions to the free energy; the bending energy tries to flatten the membrane 
and so move it away from the surface to heights where the effect of the geometrical 
heterogeneity is washed out. On the other hand, the potential energy, which, due 
to the larger surface area of the substrate, is more attractive than in the planar 
case, acts to bring the membrane in closer. As the elastic or bending energy of a 
membrane contains a rigidity term it always dominates over the potential energy 
(for a membrane exhibiting a sharp kink configuration this rigidity contribution 
even diverges - see Appendix) and thus the adhesion energy is reduced. For a flat 
substrate, the bending energy contribution is zero and this rivalry does not exist. 

This inverse dependence of U on the roughness of the substrate is perhaps 
best illustrated in Fig. |TD|, where by changing /i we move through different sur- 
face configurations whose roughness function (the variance of z s ) passes through 
a minimum. The adhesion energy can be clearly seen to mimic such behavior. In 
Fig. f|, the importance of the bending energy is again shown; with U decreasing 
as the rigidity, k, increases. 

Figures || and [11] describe the extent to which the membrane penetrates into 
an indentation in the substrate surface. Generally, the membrane extends more 
extensively into a wider trench; the bending energy of the resulting configuration 
is less — the membrane is flatter, and the potential energy gain is also greater as 
the membrane is more able to straighten out at the trench base. 

As mentioned in the introduction, supported membranes are invaluable for 
the construction of biosensors. The proteins that endow the membrane with this 
properties can, unfortunately, disturb it from a favorable planar configuration. 
One way to surmount this problem is to indent the substrate with pockets which 
then act as "docking" pods for the proteins. However, this technique can backfire if 
the membrane itself changes its configuration significantly in response to the now 
non-planar surface. Our theoretical work predicts that the narrower the trench 
or pocket, Fig. |], and the more widely spaced apart they are, Fig. [n], the less 



likely such an event occurs. We also find that, within the Deryagin framework, 
increasing the depth of the trenches encourages the membrane to penetrate further 
into them. 

Finally, we looked at randomly rough or self-affine surfaces where the gen- 
eral principles outlined above still hold. Rougher surfaces (those with a larger (5 



exponent) have smaller adhesion energies, see Fig. [12 



In an accompanying paper Ul9| , we will explore non-local methods, which en- 
able progress beyond the Deryagin approximation, and also will perform some 
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exact numerical calculations. We intend to extend the method to include chemical 
heterogeneity. This then opens up the possibility of modeling chemically and ge- 
ometrically structured substrates which are becoming crucial to recent pioneering 
biotechnological research. 
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Appendix: Kink configurations in membranes and 
interfaces 

In this appendix, we consider the elastic energy cost of a membrane adopting a 
configuration containing sharp kinks. Let us, for the purpose of illustration, con- 
sider a one-dimensional membrane embedded in a two-dimensional space having 
the profile 

h(x) = ctanhT^rJ (5.1) 

with the two lengths L and c obeying L ^> c. In the limit of e — > 0, h{x) develops 
two sharp kinks 

\\mh{x) = c[29(x) - 1] (5.2) 

and becomes a step function. We can examine the effect this has on the asymptotic 
behaviors of the two adhesion energy densities (per unit area) 

Ua = {z L dx ^ 1 + h ' 2 ( 5 - 3 ) 

and 

"* = ^£ dl (TT^)T (5 ' 4) 

Equations ( |5.3|) and ( |5.4|) provide the contributions to the free energy from the ten- 
sion (interface-like) and the rigidity (membrane-like) terms in (|2.2j ), respectively. 
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The interfacial elastic energy cost is simply proportional to the length of the 
membrane in the step configuration 



Ur, 



a 

2L 
~ <j 



2{L + c) 



(5.5) 



where we have taken the limit of e — > 0. This is clearly finite. 

To find the asymptotic limit of the rigidity contribution, we notice that 



2c 2 eK fL tanh 2 cosh 4 ( 
u K = — tt— I ax 



L 5 



+ (|) cosh" 4 (f L 
is an even function and so, writing t = tanh (^H, we have 



e 3 L 4 Jo 

which can be integrated by parts 



dt 



t 2 (l-t 2 



(5.6) 



(5.7) 



2k 
We 



1 



tanh(l/e) 



dt 



c_ ( l-t' 2 
L 



As 



ri- t anh*(lA)l 



the integral in ( |5.8| ) is well-defined and finite. 
Therefore, for small e, we have 



Ur. 



a 



and 



(5- 



(5.9) 



(5.10) 



and so the tension or interface-like contribution is constant while that from the 
rigidity (membrane-like) diverges - an interface (k = and o ^ 0) can have a 
configuration with a sharp kink while a membrane cannot. 
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K 


= 35T 


a = 1.7 x 10~ 5 Jm~ 2 


J = 38 A 


P = 





n 


~ 8.07 x 10 7 m- 1 


b ~ 5.47 x 10 4 


a- 1 = 2.2 A 


T = 


: 4.1 x 10~ 21 J 


WEAK ADHESION 


A 


= 8.67 x 10~ 22 J 


a ~ 28.5 A 


fto - 11.85a ~ 338 A 


U : 


- 1.1 X 10~ 4 cr 




- 32.25a 


v ~ 8.04 x 10" 6 a- 2 o- 


£ CT ~ 352.62a 


& = 


t 106.64a 


X 


~ 0.302 


r?+ ~ 0.0309a- 1 


~ 0.0028a- 1 


L ~ 


259.1a 






STRONG ADHESION 


(supported membrane) 






A 


= 2.6 x 10~ 21 J 


a ~ 49.3 A 


h ~ 0.61a ~ 30 A 




- 0.298cr 




- 18.62o 


v ~ 22.85a- 2 cr 


£ CT ~ 0.21a 




t 1.97a 


X 


~ 9.44 


7]+ ~ (0.359 + 0.357i)a- 1 


t?_ ~ (0.359 - 0.357i)a- 1 




22.3a 



Table 1: The various parameters, chosen and calculated, for the planar system 
and for both weak and strong adhesion. For definitions see text. 
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Figure 1: A membrane adhering to a rough surface. The reference p-plane is 
shown as a dashed line. The height of the lower membrane lipid leaflet and the 
surface, measured from this plane, are denoted by h(p) and z s (p), respectively. 
The membrane thickness is 5. 
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Figure 2: A plot of the various interactions described in Sec. 0, with parameter 
values given in Table [l], f° r a membrane weakly (a) and strongly (b) adhering 
to a planar substrate. Here all potentials are measured in units of the tension, 
a = 1.7 x 10~ 5 Jm -2 , and lengths in terms of a. The system can be seen to 
equilibrate at ho ~ 12a and ho ~ 0.6a for weak and strong adhesion, respectively. 
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Figure 3: A plot of the adhesion energy versus the rescaled wavenumber aq of a 
sinusoidally corrugated surface. For weak adhesion (a), the surface amplitude is 
set to c = 5a ~ 143 A while for a strongly adhering membrane (b) it is much 
smaller, c = 0.16a ~ 8 A. All other parameters are given in Table [I]. For these 
values of c, U soon reaches a small (< 1) asymptotic value which implies that our 
perturbation theory remains valid for all q. 
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Figure 4: The adhesion energy for a weakly (a) and strongly (b) adhering mem- 
brane above a corrugated surface. This has c = 5a ~ 143A and q = 0.02/a for the 
former and c = 0.16a ~ 8A and q = 1/a for the latter. Increasing £/a corresponds 
to increasing k/A. All other parameter values are given in Table 0. 



33 




Figure 5: The membrane profile h(x) for strong adhesion to a solid surface (thick 
line) broken by a single trench of width d = 3a ~ 148A and depth Xd with A = 0.2. 
The inset region shows a blow up of the profile with scales chosen so as to emphasize 
the overshoot present. The local adhesion energy U/Uq ~ 0.51. All lengths are 
stated in units of a without explicit statements for clarity. Parameter values are 
given in Table [I| (strong adhesion) and only the lower lipid leaflet is shown. 
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Figure 6: The height of the membrane 5h(0) = h(0) — ho at the center of an isolated 
trench of width d and depth Xd, where A = 0.2. See Table [I] (strong adhesion) for 
choices of the other parameters. 
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Figure 7: The relative adhesion energy, U/Uq plotted for a trench of width d and 
depth Xd, with A = 0.2. For an isolated, non-planar perturbation of the substrate, 
a local definition of the adhesion energy needs to be used and involves the cut-off 
L (see text). All other parameters used here are given in Table [3] (strong adhesion 
case) . 



36 



12.5 




Figure 8: A typical membrane configuration h(p) (of the lower lipid leaflet) weakly 
adhering over a cylindrically symmetric pit of depth Ar and radius r. Here A = 
0.012 and r = 350a ~ 0.997 um (for further choices see Table [I] - weak adhesion). 
The local adhesion energy (with L ~ 259.1a) is U/Uq ~ 0.50. The inset shows 
that no overshoot is present as rj± are now real. All lengths are given in units of 
a. 
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Figure 9: The profile of a membrane h(x) supported (strong adhesion) above 
a periodically structured substrate with trenches of infinite length, width d and 
depth \d. We have set A = 0.22, \i — 10 and d = 2a ~ 99 A in the above. The 
membrane configuration (thin line) follows that of the substrate (thick line) but at 
a much reduced amplitude (due to the effects of rigidity and tension). All lengths 
are shown in units of a. 
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Figure 10: For A = 0.12 and d = 2a ~ 99A, the relative adhesion energy U/Uq 
in the strong adhesion case is plotted against /x, the periodicity parameter. For 
the two flat limiting surfaces, fi = 1 and \i = oo, the decrement AU is zero and 



U — Uq. The minimum in the curve occurs near \i = 2 as described in Sec. 4Jj 
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Figure 11: The height deviation of the membrane in the middle of a trench (strong 
adhesion), 5h(0) = h(0) — h , as a function of /i for fixed A = 0.12 and d = 2a ~ 
99A. 
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Figure 12: The relative adhesion energy for a weakly adhering membrane above 
two different self-affme surfaces plotted against (, the crossover length; the full 
curve is for (3 = \ while the dashed is (3 = 1. Here, 7 = 4a, that is 7 ~ 114 A and 
all other parameter values are given in Table || (weak adhesion). 
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